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A mechanical electroscope based on a change in the resonant frequency of a cantilever one micron 
in size in the presence of charge has recently been fabricated. We derive the decoherence rate of a 
charge superpositon during measurement with such a device using a master equation theory adapted 
from quantum optics. We also investigate the information produced by such a measurement, using 
a quantum trajectory approach. Such instruments could be used in mesoscopic electronic systems, 
and future solid state quantum computers, so it is useful to know how they behave when used to 
measure quantum superpositions of charge. 

PACS numbers: 



^ '. I. INTRODUCTION 

(N : 

^ ■ As devices for processing and storing information become smaller the demands on the readout technology become 
ever greater. This is especially true for proposed solid state quantum computers which store information in various 
' quantum degrees of freedom (qubits): in quantum dots 1], nuclear spin0,|3|, superconducting islands 4] and persistent 
' currents 5], to cite just a small sample. 

Kane has proposed storing a qubit in the spin of a single phosphorous nucleus implanted in silicon. In his original 
readout scheme, this was coupled by the hyperfine interaction to the spin of the donor electron bound weakly to the 
nucleus. A surface gate would then draw the electron towards an adjacent ancilla donor, to which it might tunnel, 
^ ■ producing a doubly charged state. Under appropriate bias conditions, this transfer can only occur if the nuclear 
spin of the qubit is oriented opposite the ancilla. 

A spin measurement is thus reduced to detecting the transfer of a single electron charge to the ancilla. This can be 



o 



\ done by a sensitive electroscope such as a single electron transistor However, the techniques used for fabricating 
. microelectronics have recently been adapted to build mechanical structures at micron and even nanometre scales 
^-H ' and mechanical electroscopes sensitive to small numbers of electrons have been constructed We will consider how 
, effectively such devices might perform the measurements required for quantum information processing. 
■ Classical treatments of measurement sensitivity assume that the observable being measured has a definite value, 
^ ' which influences the measuring instrument in a definite way. The only question is how much data we must gather 
^ ' to reliably distinguish this effect from other influences on the apparatus, which produce noise. Once we know the 
' . size of the effect we wish to distinguish, and the level of noise in the system, some elementary statistics tell us the 

' integration time required for a reliable measurement. 
Q ' This assumption does not hold when we measure an observable of a quantum system. If the system is in a 
^ . superposition state, the observable will not have a definite value until some sort of measurement is carried out. Any 
^ interesting quantum information device will produce such superpositions. The process by which the superposition is 
. ^ , reduced so that the observable has a certain value imposes a minimum level of noise in the measurement, which might 
' be increased by the same sources of technical noise that affect measurements of classical systems. 
] In the proposed readout scheme for the Kane computer, a donor electron is induced to tunnel between two phos- 
. - - 1 phorous nuclei, depending on the state of the nuclear spins. In general, the nuclei are in a superposition of a state 
which would permit tunneling, and one which would prevent it. After this tunneling has occured, the electon is left 
in a superpostion of two position states, each localised on one nucleus. It then interacts with the electroscope, and 
in general with other degrees of freedom in the crystal lattice, with the result that we see it become localised on one 
nucleus or the other, so that the electroscope gives a definite signal that the charge is present or absent. 

Note that we are not discussing an ensemble of quantum systems subject to a single measurement, but rather 
a single quantum system subject to a dynamical measurement process. In such a situation we need to be able to 
describe the instantaneous conditional state of the measured system as the measurement results accumulate. This 
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is quite different from the usual situation that prevails in condensed matter systems, where typically a measurement 
is made on a large number of (almost) identical constituents undergoing quantum dynamics, and the measurement 
results are already an average over an ensemble. Fortunately mathematical techniques (known as quantum trajectory 
methods) are available to describe the conditional dynamics of a single quantum system subject to measurement with 
added noise, and these methods have been applied with considerable success to experiments in quantum optics and 
ion traps0- Recently such methods have been apphed to mesoscopic electronic svstemSjlQ. UlL IT^ . 



II. THE MECHANICAL ELECTROSCOPE 



The operation of a micro- mechanical electroscope is shown schematically in Figure^ The active part is an electrode, 
mounted on a cantilever no longer than 1 /xm, which is set in motion near the charge to be measured. The electrode 
is held at constant potential, so that its motion with respect to the unknown charge induces a flow of charge between 
it and its voltage source. The induced charge gives the electroscope electric potential energy as well as elastic, and 
changes its resonant frequency. If we envision the electroscope being used to readout a qubit in a quantum computer, 
there will be two charge states we wish to distinguish. We will denote the difference between the resonant frequencies 
of the cantilever in these two states by 6uj] it is determined by geometry and the mutual capacitance between the 
electrode and the measured charge distribution. 

We will assume the mechanical motion of the cantilever is elastic and treat it as a simple harmonic oscillator. Then 
its motion, including the capacitive coupling to the target charge, is described by a harmonic oscillator Hamiltonian 

H — h{iLjo + Sco ni)c' c (1) 

where loq is the resonant frequency of the cantilever in the absence of surface charge, and c the annihilation operator 
for its oscillation. The observable ni will be defined shortly. 

During readout of a Kane computer, a single donor electron may occupy a bound state around either of two 
adjacent nuclei. We will denote these distinct spatial states by jV') and \(f>). Only one state (suppose \ip)) couples to 
the electroscope - this is how we can distinguish them. 

During readout, the surface gates will be configured to produce tunneling between the two nuclei, depending on the 
state of the nuclear spin qubits. This entangles the charge states with the qubit states |t) and ||). We will denote 
the combined states by |0) = |t) ^ and |1) = ||) (g) according to the number of electrons interacting with 
the electroscope, which we will represent by the operator ni = |1)(1|. In general the measured qubit will be in a 
superposition state, so the total state will take the form 

1*) =a|0)+&|l) (2) 

Table 12 gives numerical parameters for a cantilever electroscope fabricated in 1998. The frequency and operating 
temperature of this electroscope meant that themal noise completely dominated any quantum effects. Besides lowering 
the temperature, this could be changed by using a cantilever with a higher resonant frequency, and such devices have 
been fabricated. However, the sensitivity of the electroscope depends on the frequency changing significantly when 
change is present, and this might not be the case in higher frequency cantilevers. 

We note that the interaction Hamiltonian commutes with the number operator ni. Furthermore, in the absence 
of tunnelling, the free Hamiltonian for the charge state itself is proportional to the square of the charge (capacitive 
electrostatic energy) and itself commutes with the charge number operator. In the presence of the measurement 
the number operator is thus a constant of motion. Such a measurement is known as a quantum nondemolition 
measurement . Number eigenstates are not changed by the coupling to the apparatus, and moments of the number 
operator are constant in time. On the other hand any state that is initially a coherent superposition in the number 
basis will be reduced to a mixture diagonal in this basis, a process known as decoherence. In an ideal quantum 
nondemolition measurement, the probability distribution for observed results at the conclusion of the measurement 
should accurately reflect the intrinsic probability distributions of the quantum nondemolition variable in the quantum 
state at the start of the measurement. 

This model, where the electroscope performs a QND measurement of the coupled charge, is idealised. If such an 
electroscope was used to measure any interesting device, the motion of the cantilever would disrupt the distribution 
of charge being measured. The nature and extent of this disruption would depend on the electrical properties of 
the system being measured; for the Kane computer, determining these is an unsolved problem in atomic physics. 
In general, back action (and interference from sources unrelated to the measurement) imposes a time limit on the 
measurement, after which the charge state will have been disrupted and the results will be meaningless. The results 
of this paper determine whether the electroscope can measure the charge with the necessary precision within that 
time. 



To detect the change in resonant frequency, we must set the cantilever in motion with some driving mechanism. In 
the device described in Tabled this was supphed by driving an ahernating current through a wire on the cantilever 
in the presence of a magnetic field. The current induced by the field in another wire was used to monitor the response 
of the cantilever to the driving. 

However, the details of the driving are not important. As long as the cantilever is coupled weakly to the driving 
system and is not damped so strongly that its state changes significantly over the period of its vibration (in other 
words, it has high finesse), the effect can be described by a Hamiltonian. In the interaction picture this takes the form 
h£{c + c'), where £ is the strength of the driving in units of frequency. If the finess of the cantilever is low, noise from 
the driving system affects its motion significantly, and the dynamics due to the driving can not be approximated by 
a Hamiltonian. 

The frequency shift could be detected in a number of ways. We could sweep the driving frequency and monitor 
the amplitude of the oscillations. Or else we could drive the oscillator at a constant frequency uj, and then detect 
the change in phase of the oscillation due to the shift in resonance frequency when a small charge is coupled; this is 
the method analysed in this paper. We will assume that if the charge state is |0), the cantilever will be driven on 
resonance; if it is the change 6uj in its resonant frequency will cause its phase to differ from that of the driving 
force. The rate of change of the phase of the output current with frequency of the driving is greatest when the 
cantilever is driven near its resonant frequency. 

We will measure time by the inverse damping rate 7~^. Then, defining a dimensionless driving strength E = £/j 
and a detuning A = Suj/j, the Hamiltonian for the coherent driving, in the interaction picture, is 



In reality of course the mechanical oscillations of the cantilever will be subject to frictional damping, and accom- 
panying mechanical noise. The rate of energy dissipation is specified by the quality factor, Q, which is the ratio of 
the resonance frequency to the width of the resonance. For linear response, this gives Q = loq/jm, where 7m is the 
decay rate of energy due to mechanical dissipation. Roukes et al.flj have measured quality factors up to 2 x lO'*. 
With such quality factors and resonance frequencies approaching GHz, these devices are approaching low quality 
optical resonators. So we will treat the effect of mechanical damping with the master equation methods of quantum 
optics. These methods assume that the coupling of the resonator to the dissipative degrees of freedom is sufficiently 
weak^,^^. Specifically we assume that ^ "^o, kT/h. 

Under these assumptions the coupling between the oscillator and the thermal mechanical reservoir is[r^ 



where a{t), a^(t) are bosonic reservoir operators. The state of the reservoir will be taken to be that of a Planck 
thermal equilibrium density operator with temperature Tm- 

We now consider in more detail the mechanism by which the small changes in resonance frequency induced by 
the proximity of a target charge are transduced. This may be done 8] by fabricating a wire loop on the mechanical 
oscillator and placing the whole apparatus in a strong magnetic field. As the mechanical oscillator moves, an induced 
EMF is set up in the loop and we may measure the induced current. When the current for the driving circuit is such 
as to drive the mechanical oscillator at its resonance frequency, the induction current is out of phase with the driving 
current. However when a small target charge shifts the resonance frequency of the oscillator, the induced current 
shifts in phase with respect to the driving current. We can detect this phase shift by an electrical comparison of 
the driving current and induction current. This is essentially homodyne detection in which the driving current plays 
the role of a local oscillator. Unfortunately this electrical transduction of the mechanical motion introduces another 
source of noise for the measurement. 

The induction current is coupled into an external amplifier circuit which can be treated as a bosonic reservoir, 
with some non zero noise temperature^^, Te- The readout circuit variable coupled to the cantilever is the current 
operator i{t) in the readout circuit. We will assume that the coupling is linear in the current and coordinate degree 
of freedom of the cantilever. Under standard assumptions the interaction between the mechanical oscillator and the 
readout circuit is described by the interaction picture Hamiltonian, 



where r{t) — 6(t)e*'^''* with the actual current in the circuit given by i{t) — y^'hu!o/2Lzo{b{t) + b^{t)) , L being the 
inductance per unit length of the transmission line, and zq the quantisation length. We will assume that the readout 
circuit reservoir is bosonic and also in thermal equilibrium at some temperature Te ■ 

Using the interaction Hamiltonians for the reservoir coupling (Equations 0] and , we may obtain the Heisenberg 
equations of motion for the oscillator and reservoir variables. Using standa7d techniques flTS^ , the reservoir variables 
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(4) 



(5) 



may be eliminated to give a quantum Langevin stochastic differential equation describing the dynamics of the oscillator 
amplitude 



^ = -iSuja- i£ -^a-^a + ^/jm ainit) + ^/^hn{t) (6) 

where ain(t), 6in(i) are the quantum noise sources for the mechanical and electrical reservoirs respectively. These 
noise terms are defined by correlation functions, which are Fourier transforms of 

(ain(t)) = (feinW) =0 (7) 

(a^HainH') = n{LO,TM)S{LU-Lo') (8) 

(ainM4,H') = in{L0,TM) + l)6iLu-Lu') (9) 

(^t^H^inH') - n{u;,TE)6{uj^uj') (10) 

(^bMbli^y) = {h{u;,TE) + l)S{u;~u;') (11) 

where 

n{uj, T)^^ (coth {hLu/2kBT) - 1) (12) 

Note the equation explicitly includes a friction term (proportional to 7b) that arises form the electrical coupling to 
the readout circuit. The steady state average amplitude a„ = (a(t))t_^ooj is given by 

(13) 



(tm + Je) + '2iSu;n 



The actual measured quantity is the current in the readout circuit, that is to say the readout variable is an electrical 
bath variable, 6out at the output from the system interaction. The output amplitudes for both the mechanical and 
electrical baths are related to the input variables for these two baths and the amplitude of the mechanical oscillator 

byia 

aout(i) = ^/lM a{t) ~ ain{t) (14) 

bontit) = -iy^ ait) - bin{t) (15) 

The average value of the electrical readout amplitude in the steady state is then found using equations Equation El 
and EauationllSI 

(bout) = Vte a« (16) 

where q;„ is given in Equation 1131 We see that the steady state amplitude of the cantilever, and hence the output 
electrical signal undergoes a change in phase and amplitude, see Figure |5| If we monitor the component in the 
imaginary direction (that is, in quadrature with the driving signal, E) we will have maximum sensitivity to this 
change in phase. Furthermore it is desirable to have E as large as possible so that small changes in phase translate 
into large changes in the quadrature. 

We can now proceed to calculating the noise pow er sp ectrum for the measured current. The calculation is analogous 
to that for a double-sided cavity given in reference|l3. We now do not work in the rotating frame but return to the 
laboratory frame. The Fourier component of the output operator for the current is given by 

[ C^''^^'^''' ) - i{ujo - u) - i6ujni]bin{uj) ~ iy/^8{uj) + ^■^ElMttiniuj) 
[{ ^ 2 ) + i{i^o - + i Suj ni\ 

where £{uj) is the Fourier component of the driving amplitude. If the driving is noiseless and monochromatic, 
£{io) = £5{ijj — ujd)- However in reality there would be some noise in the driving amplitude derived from the electrical 
noise in the driving circuit. We will treat this as entirely classical. 

Equations 1131 and 1161 suggest that the signal will appear in the quadrature of the current out of phase with the 
driving force, defined by 

X2,out{t)^i{blAt)-b,^t{t)) (18) 



with Fourier components X2^outi^)- The measured power spectrum is then given by the correlation function, 

S'2,o«t(t^,t^') = {X2,outii^), X2,outiuj')) (19) 

Using the specified states for the electronic and mechanical noise operators, we find, 

S2^out{io,Lo') - [\B{Lo)\^{2n{LU,TE) + 1) + \AiLu)\^{2n{LU,TM) + 1)] S{lu ~ lu') (20) 

where 



B{co) 



2^±^+i((c^-c^o)+ Stun,) 



1£^±21M. + i ((w - uuq) + 8uj ni) 



To estimate the signal to noise ratio (SNR) we evaluate the spectrum at the driving frequency (that is to say, at the 
central Fourier component of the coherent driving); 

"(701)2 ^ (^^^^)2l (2n{^^TE) + 1) + 7i37M(2n(w,rM) + 1) 
5(^o) = : , n (21) 



Eg uat ions i 1 31 1 and 1 1 8lshow that the magnitude of the Fourier component of the mean signal at the driving frequency 
is given by 



The signal is a sharp peak at w = = wq, in which there is a noise power S{ijJa) per root Hertz. So the SNR per 
root Hertz is |(^2,o«t(wD))P/S'(a;o), or 

[(7M + IeY + 4 <5u;2 [(^^ _ ^^)2 + 4 5^2 (2fj(^, t^) + 1) + ^ElM{2n{uj, Tm) + 1) ^ ^ 

If the SNR required for the measurement is SNR^, then we must average over noise for a time t such that SNR^ — 
SNR/\/t. If we set ni = 1, so we are measuring the charge on one electron, the sensitivity is then e^/t = e SNR^/SNR. 

III. UNCONDITIONAL DESCRIPTION OF THE MEASUREMENT. 

When we measure a quantum system, we bring an extremely large set of independent observables of our instrument 
and its environment into correlation with the measured system observable. The environment of the electroscope 
has two distinct components. Firstly there is the environment associated with the mechanical oscillator, which is 
responsible for mechanical damping and noise. Secondly there is the environment associated with the electrical 
readout, which is responsible for Johnson-Nyquist noise in the electrical circuit, and ultimately provides the measured 
result. However, we are interested in what the measurement tells us about the system, not in the exact quantum 
state of the instrument and its environment. Useful instruments must operate independently of the detailed state of 
their environments. 

There are two ways to describe the partial state of the charge and oscillator. Firstly we can ignore the results of 
the measurement and average over states of the environment completely. In this case the evolution of the charge 
and oscillator is described by a master equation. Effectively we are averaging over the ensemble of partial states 
distinguished by different measurement records 

Secondly, we can ask for the conditional states of the charge and oscillator, given a particular measurement record. 
Each member of the ensemble of partial states is associated with a distinct measurement record of the instrument. 
For it to be an effective measurement, observers must be able to distinguish the states of the instrument. In other 
words the charge must end up correlated with some simple macroscopic quantity, like the current in a wire or the 
position of a pointer on a scale. It is then possible to ask for the particular partial state of the measured system that 
is correlated with a known pointer value. In other words we need to be able to specify the conditional state of the 
system given a readout of the instrument variable that distinguishes different charge states. This is the conditional. 



or selective, description of the measured system. Of course if we average over the readout variables, we must obtain 
the unconditional description of the system. 

We begin with the unconditional description of the measurement. The dominant sources of excess noise that limit 
the quality of the measurement are the thermal mechanical noise and thermal electrical noise on the readout circuit. 
Under certain Markoff and rotating wave assumptions |0, the explicit states of the mechanical and electrical 
reservoirs may be traced out. This leaves the following master equation for the density operator of the composite 
system of charge and cantilever, 

p(t) = -i[ii{enifc^c,p\-i[E{c + c^),p\ (24) 
+ ^ -i^{n^ + l)V[c]p + -i,n^V[c^]p (25) 



^M.E 



where the superoperator T) is defined by 



I?[c]p — cpc^ — -(c^cp + pc^ c) (26) 



This can be written in a more standard form 

p = -i[p{enifc^c,p\^i[E{c + c)),p\ (27) 
+7(n+ l)X>[c]p + 7nI?[c'l']p (28) 

where 7 = 7^ + 7B: and n = {jMn{uj, Tm) + lEn{uj, Te))/j. 

We will begin solving this master equation by separating the dynamics of the cantilever and the charge. As before, 
we assume there is only one charge in the system, and consider the charge states |0) and |1). We can decompose p 
into a 2 X 2 matrix of cantilever operators 

p = i|0)(0|+B|l)(l| + Z|0)(l| + Zt|l)(0| (29) 

Since p is Hermititan, we need only three cantilever operators. A, B and Z . We can now decompose Equation 1271 into 
three independent equations involving only cantilever operators: 

rl A 

^ = -i[E{c + ct), A] + {n + l)V[c\A + nD[cf]A (30) 
^ = -i[E{c + c^) + Ac^c, B] + (n + l)V[c]B + nD[c)]B (31) 

— = -i[E{c + ct), Z] + iAZc^c +{n+ l)V[c]Z + nV[c^Z (32) 

As before, we are now measuring time relative to the damping time I/7. 

If we measured the state of the charge by means other than the cantilever, the state of the cantilever immediately 
after the measurement would be B if the charge were present, or A if it were absent. Hence A and B must be 
density operators, and Equations 1301 and 1311 have the form of master equations for a damped harmonic oscillator. 
Such equations, and their solutions, are familiar to quantum opticians. The stable solution is a displaced thermal 
state, which can be written 

p = (1 - e-^(*)) D{a{t)) e-^W^'« D^ait)) (33) 

, where D{a) is a displacement operator exp(Q;c^ — a*c), and in the steady state A = Ucuo/kbT. In the limit of low 
temperature, kT tkoQ, this becomes a coherent state |q;)(q!|. In the steady state, the cantilever has as many thermal 
phonons as a resevoir mode with the same frequency, i.e. e^^ = n/(n + 1). Its coherent amplitude ao reaches a 
balance with the driving and damping after a time around 2/7: 



a(i) = aoe-''*/'-— (l-e-'^*/2) 

K 



(34) 



1 n=0 , . 

1 + 2i A n=l 

During measurement, the cantilever states A and B are displaced thermal states with distinct coherent amplitudes. 

As the measurement proceeds, we expect the charge state to evolve from a coherent superpostion of |0) and |1) 
to an incoherent mixture; in terms of our decomposition, we expect the off-diagonal term Z to decay with time. An 
operator of the form 

Z = z{t)D{a) exp(-Actc)L>t(/3) (36) 

where z{t) is a (possibly complex) amplitude, solves Equation 1321 if a, /3, A and z obey the following differential 
equations: 

^ (n l)f - {2n + I + + n (37) 

at 

^ = {-i/\ + {n + l)l~n-\)a-iE{l-l) (38) 

'^^{{n + l)l-n-\)h + iE(l-l) (39) 

— ^^iE{a-h) + {n + \){l + ab-l) + l (40) 
at 

Here I = exp(-A), a = a-ip,b = P* - la*, and k = logz + /a*/3- i(|ap + |/3p) 

In general, these equations can be solved numerically. However, there are some special cases where we can get 
interesting information analytically. First we consider the zero temperature limit, where the off diagonal term Z is 
a projector z|q;)(/?|. The amplitudes a and /? are the amplitudes of the diagonal terms given by Equation 1351 and 
z is a complex amplitude. Once a and (3 have reached their steady state, the trace of the off-diagonal term, decays 
exponentially with a rate |a — /?p/2. 

If we assume the detuning A is small, and hence \f3\ « |q;| ~ 2E, The difference between the steady state amplitudes 
of Equation 123 is 

Cleland and Roukes give enough information about their devices for us to calculate this explicitly P| . Using the data 
m Table n we can calculate a from the definition of the annihilation operator for a torsional pendulum 



^ — (0max) = 5.3 X 10*^ (42) 

Znujo 

The normalised detuning can be calculated from the frequency shift per electron and the measured quality factor: 

A = Suj/j = = 2.4 X 10-4 (43) 

The decoherence rate is then 3.2 x 10^ 7, or 8.1 x 10^ s^^. 

As n increases from zero, the amplitudes a and /? for the off diagonal operator Z are reduced, as shown in FigureEl 
The initial decay of z{t) is shown in Figure^ and the steady decay rate, i.e. the hmit of \z'{t)/z{t)\ when t I/7, 
in Figure|3| At low temperatures (below 130 mK), the increased thermal noise from the bath causes Z to decay more 
rapidly as the temperature of the bath is increased. Contrary to expectations, the steady decoherence rate of the 
charge superposition decreases as the bath temperature increases above 130 mK. The extra thermal noise increases 
the overlap between the oscillator states corresponding to the presence and absence of charge. 



IV. CONDITIONAL DESCRIPTION 



We now turn to the correlations between the charge and the resevoir system. These are important because we must 
be able to distinguish the results corresponding to different charges to make a measurement of the charge at all. They 
can be studied most simply using quantum trajectory theory, which associates charge states with possible observed 
states of the apparatus 20]. 

We will assume we monitor the current in the electrical resevoir; this is equivalent to an optical homodyne measure- 
ment j2ll |. The inferred state of the charge as such a measurement proceeds is governed by a Wiener process, which 
is generated by a stochastic increment dW. The average of dW over the ensemble of possible measurement results 
is zero. Since the deviation of the Wiener process represented by dW increases proportional to ^/t, the average of 
(dW)'^ is dt. The simplest way to manipulate such differentials is to modify the chain rule, to give what is know as 
Ito calculus. 

Given a particular measurement result, labelled by a Wiener increment dW, the evolution of the charge and 
cantilever is 

= ii^'' - I - Kl> ^ + (1)0 + - (f » '^^) 1^) (44) 

When we insert the charge and cantilever Hamiltonian, and normalise time by the damping rate as before, this 
becomes 

dm = [-i{E{c + ct) + Anc^c)dt - \ (^c^c - 2 (|) c + (|)' j dt+(c- (|)) dW^ ^ (45) 

When a particular function dW is selected from the Wiener ensemble, this can be solved to show the evolution of a 
pure state These states form an ensemble with density operator p. Of course p can be decomposed into many 
ensembles, so the evolution generated by Equation 1451 is not unique. The details are given in Carmichael[2(Tj . 

Mixed states of the cantilever and charge must be written in the form of Equation EHl However pure states can 
always be written as 

1^) = \A) (g) |0) + \B) (g) |1) (46) 
as before we will assume the state of the cantilever is initially coherent, so that 

m = p\aO) + q\f3l) (47) 
The differential of a scaled coherent state q{t)\P{t)) is 

d{qm - (^dq - ^qd\f3f^ \l3) + q$ dt \P) (48) 

comparison with Equation 1451 gives Equations 1341 for the evolution of a and /3 as before. Some Ito calculus manipu- 
lations show that 

d\q\' - \pq\\{x)^ ~ {x)p)dW (49) 

where (x)^ is the expectation value of the amplitude quadrature x in a coherent state |a), which is just 2Rea. The 
normalisation of Itp) requires that d\p\^ = —d\q\'^. 

We need to compare the gain in knowledge shown by this trajectory picture to the decay of coherence modelled 
by the master equation. The results of the measurement are the probabilities |pp and the pure state which the 
observer will infer from these has a density operator 

p = bno)(oi + i^nixii + b<zi(ii)(oi + io>(ii) (50) 

The off-diagonal terms in this have magnitude \pq\; we can average over dW to see the behaviour of the density 
operator for the ensemble of measurement results. 

Some more routine Ito calculus gives the evolution of this: 

d\pq\ = -\pq\ Q {{x), - {x)^)'dt + \ (bp - Igp) ((x)^ - (x) J dW^ (51) 



Since the average of dW over different measurement results is zero, on average 



d\pq\= -l{{x)^-{x)J^\pq\dt (52) 

If the difference between the charges associated with states |0) and |1) is e, then in the state p\aO) + q\pi), the 
uncertainty in the charge is given by 

1 

[{{nef) - {nefy ^ e\pq\ (53) 
From Equation 1521 this decreases exponentially as the measurement progresses, at a rate 

This differs from the square root decay of classical uncertainty as measurements are averaged over time, but exponential 
decay is what we would expect for decay of coherence 13). For the device described in 8], this is almost equal to the 
decoherence rate. In real devices, thermal noise will cause the trajectory states to be mixed, however the evolution of 
such mixed states is much harder to calculate. 



V. DISCUSSION 



To estimate the time required for our measurement, we have calculated how long it takes for an initially pure 
superposition of charge states to be reduced to a mixture, and how long (in some sense) it takes us to find out which 
charge eigenstate we have been left with. While these questions are interesting in their own right (they composed 
the deepest mystery of physics for the best part of a century), it could be argued that they don't reflect the way 
measurements would be used in a real computer. 

The most that we could do with measurements on pure states is state preparation. In a coherent quantum computer 
this would be rather pointless though, since if we know the initial state we could just rotate it into the eigenstate we 
want. We carry out measurements to find out something wc don't know: in other words we apply them to mixed 
states, with a view to finding out which of the possibilities is real. 

Information theory provides tools to quantify this, such as conditional entropy and mutual information. Unfortu- 
nately calculating any of these requires knowledge of the ensemble of trajectories generated by each component of the 
mixture, and the overlaps between them. In general it is hard to find the probability distribution of trajectories; we 
usually just calculate averages. It might be worth doing this numerically, however. 

There is a more straightforward limitation to our analysis: in present day devices the thermal effects that we have 
neglected in the trajectory treatment utterly dominate the vacuum noise we have considered. Hence the measurement 
time will be limited by the need to average classical fluctuations. It is possible that future devices operating at 
higher frequencies will reduce the level of thermal noise so that quantum effects will be important. This presents the 
remarkable prospect of a solid cantilever with position and momentum known to the limit allowed by the uncertainty 
principle. 
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Operating temperature T 4.2 K /cbT = 3.6 x 10"'* eV 

Resonant frequency (jJo/27r 2.6 MHz huj^ = 1.1 x 10~* eV 

Torsional spring constant k 1.1 x 10"^" Nm 

Amplitude 0max 30 mrad 

Frequency shift per electron 5v 0.1 Hz 

Quality factor tJo/7 6.5 x 10^ 

TABLE I: Data for an electroscope fabricated by Cleland and Roukes |i 



Signal 





Electrode 



V / Charge 



Induced 
"^^ ' Charge 




FIG. 1: Operation of the mechanical electroscope. A charge trapped near the surface of some material is coupled to a cantilever 
suspended above the surface, as explained in the text. The cantilever is driven at a rate E and damped by a combination of 
mechanical friction and reaction from the electronic readout loop at a rate 7. If an excess charge is present on the surface, the 
frequency of the pendulum is increased by 5. For simplicity, the figure shows a simple pendulum, but in practice the cantilever 
would be a torsional pendulum, oscillating due to strain in the material. 



Amplitude Quadrature 10 ^Xi 
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FIG. 2: The off-diagonal element Z of the density operator is a thermal state displaced by amplitudes a and /3, which depend 
on the temperature (see Eauation l36ll . When the cantilever is coupled to a hot bath, these coherent amplitudes decrease, and Z 
approaches a purely thermal state. The values these amplitudes take in the Cleland and Roukes electroscope at temperatures 
from absolute zero up to 10 K are plotted in the complex plane, in units of the ground state fluctuations. The amplitudes of 
the diagonal elements A and B do not vary with temperature, but remain at the K values. 
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FIG. 3: When the measurement has been running for a time around I/7, and the cantilever amplitudes have reached their 
steady state, any remaining coherence between the two charge states decays exponentially. Here the rate of this decay is plotted 
as a function of temperature, for the device described in Table Q] The maximum decay rate of —3.2 x 10^7 occurs at 130 mK. 
Beyond this point the decay rate decreases with temperature, possibly because the increased thermal noise makes the coherent 
amplitude of the cantilever harder to distinguish. 
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FIG. 4: The state of the cantilever takes some time to become entangled with the charge after they begin to interact, as the 
cantilever state moves towards its steady value. After this the charge state decoheres rapidly. Here the coherence between the 

two charge states is plotted as a function of time for an array of temperatures. The cantilever is initially in a thermal state at 
the appropriate temperature. Note the the charge state has decohered long before the cantilever reaches its steady amplitude, 
which occurs after a time I/7. 



